Electron energy relaxation in the presence of magnetic impurities 
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We study inelastic electron-electron scattering mediated by the exchange interaction of electrons 
with magnetic impurities, and find the kernel of the corresponding two-particle collision integral. In 
a wide region of parameters, the kernel K is proportional to the inverse square of the transferred 
energy, K oc J 4 /E 2 . The exchange constant J is renormalized due to the Kondo effect, yielding 
an additional weak dependence of K on the energies of the colliding electrons. At small energy 
transfers, the 1/E 2 divergence is cut off; the cut-off energy is determined by the dynamics of the 
impurity spins. The obtained results may provide a quantitative explanation of the experiments of 
Pothier et al. [Phys. Rev. Lett. 79, 3490 (1997)] on anomalously strong energy relaxation in short 
metallic wires. 
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The effect of magnetic impurities on the electron prop- 
erties of a metal is drastically different from that of 
"usual" defects which just violate the translational in- 
variance of the crystalline lattice. The reason for the dif- 
ference is that a magnetic impurity brings an additional 
degree of freedom - its spin. If there were no exchange 
interaction with the itinerant electrons, the ground state 
of the system would be degenerate with respect to the 
orientation of impurity spins. Weak exchange interac- 
tion allows an itinerant electron to flip its spin in the 
course of scattering on a magnetic impurity. Such scat- 
tering, accounted for even in the lowest-order (Born) ap- 
proximation, yields an important effect of dephasing of 
the electron state. Finite dephasing time, in turn, sup- 
presses the interference corrections to the conductivity, 
thus suppressing the weak localization effect Q . 

The higher-order terms in the perturbation theory se- 
ries for the scattering amplitude reveal one more im- 
portant phenomenon. It turns out that the amplitude 
of scattering caused by the exchange interaction in- 
creases with lowering the temperature, as opposed to the 
temperature-independent scattering on a usual impurity. 
This increase is responsible for the non-monotonic tem- 
perature dependence of the resistivity of a metal, the 
phenomenon called the Kondo effect ||. 

Spin exchange between an electron and a magnetic im- 
purity may occur in an act of elastic scattering. Account- 
ing for these spin-exchange elastic processes is sufficient 
for understanding the dephasing phenomenon Q and the 
Kondo effect J|]. However, such processes do not lead 
to any energy relaxation of electrons. In this paper we 
demonstrate that magnetic impurities may also mediate 
energy transfer between electrons. If the energy transfer 
E is larger than the Kondo temperature Tk, then the 
energy relaxation occurs predominantly in two-electron 
collisions. We derive the kernel K of the corresponding 
collision integral in the kinetic equation for the distribu- 
tion function. This kernel depends strongly on the trans- 
ferred energy, K oc J 4 /E 2 . The dependence of K on the 
energies e, of the colliding electrons (measured from the 



Fermi level) is relatively weak as long as ^S> Tk- This 
dependence comes from the logarithmic in | er^ j renormal- 
ization of the exchange integral J, known from the theory 
of Kondo effect At small energy transfers, the l/E 2 
divergence of the kernel is cut off; the cut-off energy is 
determined by the dynamics of the impurity spins, which 
results from their interaction with the Fermi sea. 

The motivation for our study comes from the experi- 
ment where the relaxation of the electron energy dis- 
tribution function in mesoscopic wires was investigated. 
It was found that the empirical relation K oc 1 /E 2 holds 
in a substantial interval of energies E for Au and Cu 
wires. The data of Ref. {|,H was accurate enough to rule 
out the direct Coulomb interaction || , which would yield 
K{E) oc 1/E 3 / 2 , as a source of relaxation. 

We describe the metal with magnetic impurities by 
means of the exchange Hamiltonian: 
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where index I labels the magnetic impurities, S/ is the 
spin operator of the Z-th impurity, = S{S + 1), and r; 
is its coordinate. Free electron states are labelled by the 
wave vector k and spin index a. The Pauli matrices are 
denoted by a = (a x ,a v ,a z ). 

If the concentration n of the impurities is low enough, 
they can be considered independently. Therefore we will 
perform our calculations for a single impurity, omitting 
the impurity index I, and then will multiply the resulting 
expressions for the scattering rate by n. In this one- 
impurity problem, there is interaction only in s chan- 
nel, so we will label the participating electron states with 
scalar index k. 

In the framework of the exchange Hamiltonian (|lj) , the 
lowest non- vanishing order of the perturbation theory se- 
ries in the exchange constant J for the inelastic scattering 
amplitude is the second order: 
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In the diagrammatic representation, the amplitude is 
the sum of the diagram shown in Fig. [l] and the three 
other diagrams that can be obtained from the diagram 
of Fig. [l] by the transposition of indices 1 2 and/or 
3 «-* 4. Note that there is no summation over the initial 
or final spin states of the impurity in Eq. (Q). Therefore, 
the spin lines are not closed, i.e. contrary to Ref. 0] this 
scattering amplitude cannot be represented in the form 
of an effective four-electron vortex. 
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FIG. 1. A characteristic diagram for the amplitude of in- 
elastic electron-electron scattering mediated by the exchange 
interaction of electrons with a magnetic impurity, in notation 
of Ref. The solid lines denote electron states, the dashed 
lines denote the localized spin state. 

Performing the time integrations in Eq. (|^) , we obtain 
the standard expression for the second-order term of the 
expansion of T-matrix, 



A(kiai, k 2 a 2 , S fc 3 CT 3 , A; 4 c r 4, S') = 

(k 3 a 3 , k 4 a 4 , S'\V- — ^-V\hai, k 2 a 2l S) 
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The denominator in Eq. (||) is the energy of the inter- 
mediate virtual state, which equals ±(£fci — £fc 3 ) for two 
of the four possible pairings of the electron creation- 
annihilation operators in Eqs. (^) and (||) (one of these 
pairings is shown on Fig. |l|), or ±(^ 1 — £fc 4 ) for the other 
two pairings. The spin structure of the scattering ampli- 
tude can easily be found from Eq. <^). In a scattering 
event, spins of one or both participating electrons must 
flip, with the corresponding change of the impurity spin. 
In this paper we are interested only in the relaxation of 
the electron energy distribution, and assume that the sys- 
tem does not have any spin polarization. Therefore we 
need to calculate only the total cross-section of scattering 
into all possible spin states, averaged over the initial spin 
states of the impurity and two electrons. After this aver- 
aging, the terms proportional to -£k 3 )(£ki -CfcJ] -1 
drop out. Finally, we get the rate of scattering of two 
electrons with energies e± and e 2 into states with ener- 
gies £3 and £4: 



r(£ 1 ,£ 2 ;£ 3 ,£4) = --S{S+l){Jvf 

4 v 



(4) 



(£1 - £3) 2 (ei - £4) 



S(si + £ 2 - £3 - £4) 



Substituting this expression into the collision integral Q 
and taking into account the symmetry e\ <-> e 2 , £3 *-* £ 4 
of the right-hand side of Eq. ([!]), we obtain the collision 
integral in the form: 

/(£,{/}) = J J dEde'K(E) (5) 

x {/(£)/(£') [l-f(e-E)} [l-f(e' + E)} 
- f(s-E)f(s> + £)[!- /( e )][l -/(£')]} , 



where the kernel 



K(E) = ^S(S 
2 v 
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depends only on the energy E transferred in the collision, 
and /(e) is the electron energy distribution function. 

The above derivation of the inelastic amplitude 
[Eq. (^)] was performed in the lowest non-vanishing or- 
der of the perturbation theory. It is known from the 
theory of the Kondo effect, that for the elastic scatter- 
ing amplitude, the calculation in the lowest order may 
be insufficient. The higher in Jv orders yield contribu- 
tions to the elastic scattering which logarithmically di- 
verge at low energies For elastic scattering, the 
leading terms in all orders can be summed up with the 
help of the renormalization-group technique |l0| . In this 
technique, the bare exchange constant is replaced by the 
renormalized one, 



J(e) = [vU{\e\/T K )Y 



(7) 



and the scattering amplitude at energy e is to be cal- 
culated within the Born approximation in J{e). The 
Kondo temperature Tr is related to the parameters of 
the Hamiltonian (fy) by 



(8) 



Tk = M"v JvD exp (—1/ Jv) 



Here D is the high-energy cut-off, and /1 ~ 1 (for detailed 
discussion of these parameters, see Ref. (|Tl|). 

Similar to the theory of elastic scattering, the lowest- 
order result (||) is valid only while hi(\ei\/Tx) ^ 1- At 
smaller energies, the two vertices of Fig. [j] acquire cor- 
rections of the form Jv[Jv\n(D /\ei\)] m in the (m+ f)-st 
order of the perturbation theory. The two vertices are 
renormalized independently from each other [ jl2| , and the 
corresponding diagrams have the same structure as in 
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Ref. 0. This approximation is justified as long as the 
energies Si of all incoming and outgoing electrons satisfy 
the condition 



]n(\e t \/T K )>l 



(9) 



The resulting non-perturbative expression for the kernel 
in the collision integral reads 



K(E;e,e') = 8n-S(S + l) 
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We would like to emphasize that the dependence of the 
kernel K on the energies of participating electrons re- 
mains weak (logarithmic) , as opposed to the strong l/E 2 
dependence on the transferred energy E. The logarith- 
mic dependence in Eq. (]l(]) is meaningful as long as the 
electron energies £j exceed the smearing of the Fermi level 
caused by temperature or a non-equilibrium electron dis- 
tribution. In the opposite case, the arguments of all the 
logarithms in Eq. ( |l0| ) must be replaced by e* /Tk, where 
energy e* characterizes the smearing. It is important to 
note that e* does not cut off the singularity in the trans- 
ferred energy E. 

The low-energy divergence of the scattering rate (Q) 
stems from the degeneracy of spin states of the impurity. 
Due to this degeneracy, the energy of the intermediate 
virtual state [the denominator in Eq. (|3|)] approaches zero 
at E — > 0. An additional condition for the divergency is 
the time independence of the average (S'\S : '(i)S k (t')\S) 
in the approximation of Eq. (||) (here j, k — x,y, z). In 
fact, exchange interaction between the itinerant electrons 
and impurity may flip its spin. The resulting impurity 
spin correlation function decays, restricting the lifetime 
of the intermediate virtual state. The corresponding de- 
cay rate cuts off the E = singularity of the kernel (10). 
The manner of decay depends on the electron energy dis- 
tribution /(e). We will discuss first the cut-off in the case 
of weak deviations of /(e) from the thermal equilibrium. 

Let us first consider the low-energy cut-off for K(E) 
at high temperatures, T 3> Tk- Scattering of electrons 
off the spin results in exponential decay of the correlation 
function, (S'\Si(t)S k (t')\S) oc exp(-\t-t'\/T T ). The im- 
purity spin correlation time tt can be evaluated with the 
help of the Fermi golden rule. Since the deviation from 
the thermal equilibrium is weak, we can replace /(e) with 
the Fermi distribution function uf{s), 



— = ^S(S + \){Jv) 2 / den F (e)[l - n F (e)\ 
tt A J 



= yS(S+l)W 2 T 



(11) 



When T is lowered towards Tk, the exchange constant in 
Eq. ( pd| ) is renormalized according to Eq. ^ with e = T. 
The resulting expression for the spin-flip rate reads: 



- = ^-S(S + l)[MT/T K )Y 2 T . (12) 
tt o 



The energy scale Ti/tt sets the limit of applicability of 
Eq. (||) and cuts off the singularity in the kernel (|lo| ) at 
E ~ Ti/tt- Note that within the limits of applicabil- 
ity of Eq. (12), the spin-flip rate satisfies the condition 



T > h/ tt > T K - 

At T < Tk, we can neglect the effect of a finite temper- 
ature on the scattering of electrons with \ei\ ^> Tk- The 
behavior of K(E) can be related to the time dependence 
of the zero-temperature self-correlation function of the 
impurity spin. At time difference \t — t'\ 3> Tl/Tk this 
correlation function decays as (GS\S j (t)S k (t')\GS) ~ 
5jk/[TK\t — 1'\] 2 , and therefore the singularity of the ker- 
nel K(E) is cut off at E ~ T K (here \GS) is the ground 
state of the Kondo problem) . 

At this point, we should mention that at energy trans- 
fers E ~ Tk, the processes with participation of three 
or more electrons must also be taken into account along 
with the two-particle scattering. The consideration of 
these multi-particle processes is an arduous task lying 
beyond the scope of this paper. Here we address mostly 
the electron energy relaxation on large Tk) energy 
scales. It allows us to limit our consideration to the two- 
particle processes accounted for by Eq. (^) , and dispense 
with the scattering events involving more particles. 

At very small energies (|£j|,T <§C Tk), however, the 
Fermi-liquid description of electrons is again a valid tool. 
The behavior of the system is described in this case by 
the quadratic fixed-point Hamiltonian, in which the four- 
fermion interaction is a least-irrelevant term ^3[|lJ] . The 
calculation of the inelastic scattering rate is then straight- 
forward, the resulting collision-integral kernel is given by 



K(E) 



1 n 



(13) 



When T = 0, the corresponding rate of inelastic electron 
scattering is h/r in = dEK(E)E cx (e/T K ) 2 - At e -> 0, 
it decreases faster than e, as it is supposed to be in the 
Fermi-liquid picture. 

Relaxation of the electron energy distribution was in- 
vestigated experimentally in metallic wires of Cu and 
Au in Refs. ^,|). In these experiments, a finite bias 
V = 50 — 500 [xeV was applied to the ends of a wire. It 
was found that starting from fairly small wire lengths, the 
electron distribution is smeared over the range of energies 
eV, instead of having two distinct steps created by the 
bias applied to the wire ends. The observed electron en- 
ergy relaxation was attributed ^|,|5| to two-electron colli- 
sions. The collision-integral kernel for E < eV extracted 
from the experiments has the form K(E) = Tl/(tqE 2 ) : 
with a cut-off at some low energy, which scales linearly 
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with eV Jig]. The value of the parameter tq was 0.5 — 1.0 



ns for Cu wires, and 0.1 ns for Au wires |16|. 

Now we discuss the possibility of such relaxation due to 
the electron scattering on magnetic impurities in wires. 
In the experimental setup, the electron distribution is 
smeared, and the width of smearing eV exceeds both T 
and the energies \ei\ of the colliding electrons. Assuming 
also that eV ^> Tk, we can simplify the kernel (flQ): 



l)[\n(eV/T K )}-^ . (14) 



The l/E 2 dependence in Eq. (|l|) persists down to 
the cut-off, which is determined by the spin-flip rate 
1 / T eV ■ The derivation of the spin-flip rate in these non- 
equilibrium conditions follows the same path which led 
to Eq. (p^), and results in: 



TeV 



~/S(S + l)[\n{eV/T K )]- 2 eV. 



(15) 



Here the numerical constant 7 ~ 1 depends on the de- 
tailed shape of the non-equilibrium electron distribution. 

Now we compare theoretical results ([l4]) and ( |l5| ) with 
the experimental observations for Au wires ||. Proper- 
ties of these samples are compatible with the presence 
of iron impurities with a concentration up to few tens 
of ppm |15|. We take the density of states in Au at the 
Fermi level v « 0.25 (eV site)" 1 and T K w 0.3K for 
Fe impurities in Au |2| . The typical value of voltage in 
the experiments Q was V ~ 0.3 meV fill . Substituting 



these parameters into Eq. (14), we see that a relatively 
small concentration n ~ 10 ppm is sufficient to repro- 
duce the experimentally measured value To ~ 0.1 ns p8|| . 
The spin-flip rate ([l5]) is the low-energy cut-off for the 
l/E 2 dependence of the kernel. This cut-off is roughly 
proportional to the applied voltage, in agreement with 
experimental observations p5| . We must note, however, 
that the lower voltages used in experiment || are close 
to the Kondo temperature, so the leading-logarithmic ap- 
proximation ||s|,|lC)|l, used in derivation of Eqs. (|l^)-(|T5|), 
may be insufficient. 

The inelastic electron scattering off magnetic impuri- 
ties must be sensitive to an external magnetic field po- 
larizing the spins in the system. The Zeeman splitting 
gfiH of the spin states prevents the impurity spins from 
changing, thus suppressing the inelastic scattering pro- 
cesses with energy transfers E < g/J.H. Measurements in 
a sufficiently strong magnetic field may elucidate the role 
of magnetic impurities in the electron energy relaxation. 

An important feature of the electron-electron inter- 
action mediated by magnetic impurities is that it is 
not translationally invariant. This is why the introduc- 
tion of non-magnetic impurities, which affects drastically 
Coulomb interaction of electrons j(| , produces only small 
corrections to the interaction induced by magnetic impu- 
rities. 



The above consideration was performed for the mag- 
netic impurities described by the one-channel Kondo 
model. However, the proportionally of the scattering in- 
tegral kernel to l/E 2 at E ^> Tk holds for an arbitrary 
number N of channels in the Kondo problem. If the ex- 
change constants J are the same for all channels, the 
kernel (^) acquires an additional factor N 2 . For the spe- 
cific case N = 2, the qualitative K(E) cx l/E 2 behavior 
was noticed in Ref. [|l9|. However in Ref. |19| this behav- 
ior was attributed solely to N = 2; this, in our opinion, 
is inaccurate. Our consideration also indicates that at 
any N, the l/E 2 divergence of the kernel K(E) is cut- 
off at small E, see Eq. (15). This is also in an apparent 
disagreement with Ref. [19], which states that for N = 2 
the divergence persists down to the lowest energies. 

In conclusion, we have shown that the exchange in- 
teraction of itinerant electrons with magnetic impurities 
can facilitate inelastic electron-electron scattering. We 
derived the kernel of the corresponding collision integral, 
and found its explicit dependence on the parameters of 
the system for a wide range of the energies of colliding 
electrons. This allowed us to perform a quantitative anal- 
ysis of the experimental results of Refs. We find 
that a very small density of magnetic impurities could 
lead to the anomalies in the electron energy relaxation 
observed there. 
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